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Abstract 

Zurek, Habib and Paz [W. H. Zurek, S. Habib and J. P. Paz, Phys. Rev. Lett. 
70 (1993) 1187] have characterized the set of states of maximal stability de- 
fined as the set of states having minimum entropy increase due to interaction 
with an environment, and shown that coherent states are maximal for the par- 
ticular environment model examined. To generalize these results, I consider 
entropy production within the Lindblad theory of open systems, treating envi- 
ronment effects perturbatively. I characterize the maximally predicitive states 
which emerge from several forms of effective dynamics, including decoherence 
from spatially correlated noise. Under a variety of conditions, coherent states 

emerge as the maximal states. 
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I. INTRODUCTION 



Decoherence which results from a quantum system's interaction with an environment can 
provide a mechanism for characterizing the transition from quantum to classical behavior for 
a quantum open system |I]J2], and has been an integral part of several programs addressing 
the emergence of classicality ||[§]. Zurek, Habib and Paz (ZHP) have characterized the 
effectiveness of decoherence in terms of a predictability sieve, and identified the maximally 
predictive states (defined as those with minimal entropy production) as the most classical 
B. ZHP considered an environment model consisting of an independent oscillator bath 
linearly coupled to the system of interest, which has been studied in the context of quantum 
brownian motion |J . ZHP demonstrated for the high temperature limit of the environment 
that the coherent states of a harmonic oscillator are maximally predictive, and that zero 
squeezing (corresponding to coherent states) is maximal for squeezed states considered at 
arbitrary environment temperatures. 

The purpose of this paper is to extend the results of ZHP to additional environment 
models. The effects of the environment are considered in the general framework of the 
Lindblad form for nonunitary evolution of a harmonic oscillator 0, correponding to the 
Markov limit. In Sec. [TT] I establish a systematic framework for evaluating the predictability 
of states using first order perturbation theory. In Sec. |T| I apply this framework to a family 
of Lindblad generators which are at most quadratic in position and momentum whose general 
properties have been studied in the literature f|. This family of Lindblad generators includes 
as special cases the Quantum Optical Master Equation |j] and Dekker's model for Quantum 



dissipation ||10|| . I consider models with environment spatial correlation effects [TT],[T2| in 
Sec. [TV]- I comment on these results in Sec. [V]. 
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II. PERTURB ATIVE APPROACH TO EVALUATING PREDICTABILITY 



The predictability sieve was introduced by Zurek || as a means of characterizing those 
states which are most stable when considering not only the dynamics of the system, but 
also including the effects (such as decoherence) of interaction with the systems environment. 
The set of states having minimum (linear) entropy production are the best candidates for 
corresponding to points of classical phase space. The linear entropy ^ for a particular density 
operator p is given by: 

,(t) = Tr[p(t)-p(t) 2 ] = l-Tr[p(t) 2 ]. (2.1) 

The effect of the environment on evolution is taken to be in the Markov regime, and I assume 
that there is no further explicit time dependence on the Liouville generator for the isolated 
system. The form of the evolution is then stationary, and we can write 

p(t) = e Lt [p} = J{t)[p\. (2.2) 

The generator L is assumed to be comprised of two parts, the evolution of the isolated system 
corresponding to L with J {t) = e Lot , and the effects of the environment AL. I consider 
the effects of the environment perturbatively in part because a wide variety of forms of 
evolution become tractable, and also because many of the effective evolution equations are 
derived using lowest order approximations such as the weak coupling limit fl3| . From formal 
perturbation theory, the evolution of Liouville generator can be written 

J(t) = J (t) + I J Q (t - T)ALJ{r)dr 
Jo 

= Jo(t) + f J (t ~ T)ALJ (r)dr + 0{AL 2 ) 
Jo 

= J a (t) + Ji(t) + 0(AL 2 ), (2.3) 

defining J\{t) as the first order perturbation. For an initial state p Q , the entropy at time t 
becomes 
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q(t) = l-Tr[J(t)\p„]J(t)\p ]] 
^l-Tr[J (t)[p }J (t)[p }} 

-2Tr[J (t)[p ]J 1 (t)[p ]]. (2.4) 

If the isolated system has unitary evolution generated by a hamiltonian H, then 

L \p] = ±:[H,p], 
J [-]=U(t)-U\t), 

U(t) = e^ Ht . (2.5) 

In this case 

1 - Tr[J (t)[p ]J (t)[p \] = 1 - Tr[ PoPo ] = ? (0). (2.6) 
The maximal states are those with minimal entropy production, that is, with minimized 

Aq(t) = -2Tr[J (t)[p ] Ji(t)[p ]\ 

= -2 f{^[U{t) Po U\t)U{t - r)AL[U{s)p U\s)\U\t - r)]}dr 
Jo 

= -2 f TT[p UKr)AL[U(r)p U f (r)P(r)]dT, (2.7) 
Jo 

where the cyclic property of the trace and the unitarity of U (t) have been used for the final 
simplification. 

The general form for a Lindblad generator is 

L[p] = h[H,p] + ^[Vjp,^] + [V^pVj], (2.8) 

3 

so that the contribution to environment interaction can be written as the perturbation 

AL[p] = ^H,p] + ^Et^^l + Wi,PV}]- (2-9) 
The entropy production can now be written 

Aq(t) = -2 [' Tr[p U\r)^-[AH,U(T)p U^T)}U(r)dT 
Jo in 

~\ J* Tr\p U\r) Y,MU(t) Po U\t), V]} 

+ [V V U{r) Po U ] {r)V^)U{r)]dT. (2.10) 
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The first term on the RHS of this equation is identically zero, from the cyclic property of the 
trace. The remainder can written in a simpler form using the cyclic property of the trace, 
the unitarity of U(t) and the identification: 

= (2-11) 

to yield 

A ^) = ^/ t ETr[{V/(r)V J (r),p o } 

-2V 3 {r) Po V^r)}dT. (2.12) 



If p is a pure state, then it is also a projection, with 

p = P=\1>)y,\ (2.13) 

with 

P 2 = P (2.14) 

and 

POP = P(ip\0\ijj) = P(0) (2.15) 
for an arbitrary operator O. Thus, for pure states, the entropy production is given by 

= l^E{({^{r)V s (r)) - (V,(r)}(Vl(r)})dT. (2.16) 
Minimization of this final quantity can then be used to determine the maximal states. 

III. LINDBLAD OPERATORS QUADRATIC IN POSITION AND MOMENTUM 

I will now apply the results of Sec. || to Lindblad operators which have {Vj} linear in 
position and momentum. This family of generators have been studied extensively in the 
literature M, and include as special cases the quantum optical master equation |§ and 



Dekker's phenomenological master equation 1 10(1 - In terms of Eq. ( |2.9|) , the operators are 
given by 

AH=^{x,p} 
Vj = cijp + bjX. (3.1) 



With the identifications 
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z j 



A = ElmK&*], (3.2) 



the perturbation on the system evolution becomes 



AL[P] = k [ 2 {x ' Php] ~ M ~ [X ' P]] 

+ ^([x,[p,p]] + [p, [x,p]]) 

(3-3) 



The particular choice of parameters Z) OT , -D pp , -D p? , A and /i Eq. (|3.3| ) determines the details of 
the evolution (i.e. evolution corresponding to Quantum Optical Master Equation, Dekker's 
master equation, etc.). It has been determined that only when 

A> ? = ~E- Re [ Q A*] = ° (3-4) 

Z 3 

will the system relax into a thermal equilibrium state 0, so I adopt this condition for the 
remainder of the paper. 

For the simple harmonic oscillator, the operator equations of motion are easily solved by 

x( T ) = U ] (t)xU(t) 

P 

= x cos(to>r) H sm(uT), 

mui 

p(r) = U\t) P U(t) 

= pcos(ujt) — mux sin(c<jr), (3-5) 



to yield 



Vj(r) = aj(pcos(ur) — mux sin(a;r)) 

+bj(x cosftur) H sinftur)). (3.6) 



Substituting Eq. (|3.(j|) into Eq. (|2.16p and evaluating the elementary trigonometric integrals 



over r: 



where 



1 Tfllil 

+/ 2 W(^(p 2 -(p) 2 )-^(x 2 -(x) 2 )) 

+fs(t)^({x,p})-uj(x}{p}), (3.7) 



/ 1 (t) = t2m(^ + ^%), 
/ 2 (i) = 2 m - Sin(2 " i) ' 2B " 2D » 



2cj v ft 2 (muohf h 

h{t) = _ 2 E^<^ _ (3. 8) 
c<j z ft, {moony 

Eq. ( |3.7D is the expectation of a c-number times the harmonic oscillator hamiltonian plus 
a second c-number constant after sqeezeing and translation by (x) in position and (p) in 
momentum. The state which minimizes Eq. ( |3.7| ) will be the corresponding squeezed and 
translated ground state, which is simply a coherent squeezed state ||14|| . In terms of the 
harmonic oscillator creation and annihilation operators a) and a, the squeeze operator is 
given by 

S(() = e §«* a2 -^ t2 ), (3.9) 
and Glauber's displacement operator is given by 

D(a) = e {aa ^ aa) . (3.10) 
The appropriate selection of the paramter a for the displacement operator 



2muj 

„ . . . hmuj . i , . . 
HB(a) = (— )i<rr>, (3.11) 

provides the necessary translation (in phase space). For simplicity, I take (x) = and 
(p) = for the rest of this section. The effect of squeezing on the anihillation operator can 
be written 

S*(t)aS(0 = fm + va\ (3.12) 

where 

C = se*, 
\i = cosh(s), 

v = smh(s)e ie . (3.13) 

The parameter s determines the amount of squeezing (s = corresponding to no squeezing), 
and the parameter 9 determines the orientation of the sqeeze axis. The squeezed harmonic 
oscillator hamiltonian is given by 

S f (C)(fU^ 2 )S(C) 
= cosh(2 s) (— J + 

+ sinh 2s sin ^ x 2 - — 

2 2m 

+ cosh(2s)cos(^)^{x,p}, (3.14) 

which can be used to identify the amount and direction of squeezing required to minimize 
Eq. ( |3.7| ). The easiest condition to extract is the direction of squeezing: 

cos(0) = |||. (3.15) 

Since f\{t) is linear in t while fsit) oscillates, the long time behavior is 9 = 0. The amount 
of squeezing can be determined by examining 
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(1 + cosh 2 (2s)) sin(0) = ^\ cosh(2s). (3.16) 

To analyze the long time behavior, it is useful to note that sin(0) approaches 1 and f2if)/ fiif) 
approaches (and therefore so does the RHS of Eq. (|3.16| )). Thus, the long time behavior 
of the amount of squeezing required for minimum entropy production is 

(1 + cosh 2 (2s)) = 0, (3.17) 

which implies s = 0, no squeezing. Thus, for times on the order of the dynamical time scale 
of the system (more than a few cycles), coherent states are the maximal states, just as ZHP 
found for their environment model. 

IV. CORRELATION EFFECTS IN ENVIRONMENT NOISE 

In this section I wish to consider environment models which include the effects of finite 
correlation lengths in the environment |pT|,12|. From a strictly phenomenological point of 



view, Quantum Mechanics with Spontaneous Localization can be included by virtue of the 

effective form of the dynamics, although this is actually a fundamental modification of 

Quantum Mechanics fll5| . Many of these models can be written in the form 

dp(x,x';t) TT . 

= Hamiltoman + Dissipation 

at 

+ g(x,x')p(x,x';t), (4.1) 

where the decoherence term which we shall focus upon can be expressed in terms of the 
correlations of a classical fluctuating potential V(x, t) with 



g{ x \ y) = ( c (^; x ) + c (y; v) - 2c O; y)), ( 4 - 2 ) 



where 

(V{x, t)V(y, r)) av = c(x; y)S(t - r). (4.3) 

For simplicity, I consider only a homogeneous and isotropic environment for which c(x; y) = 
c(x — y), and g(x;y) = g(x — y). I will also restrict my attention to weak dissipation 
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, and consider only the evolution due to the unperturbed hamiltonian and the (spatially 
correlated) noise term. 

The Lindblad form can be used to represent the noise term with 

{Vj} = a{k)e lkx (4.4) 

and replacing the discrete sum over j in Eq. fl2.9|) with an integral over k. The evolution in 
this case can be written 

dp ( x £'> f 2 = Hamiltonian - i j dk\a(k)\ 2 (l - e ik ^- x,) )p(x,x';t), (4.5) 

so that |a(fc)| 2 and c(r) are Fourier transform pairs: 

c (r) = 1 1 \a{k)\ 2 e ik{r) dk. (4.6) 



Thus a noise term with short correlation length scales will have a narrow c(r) and a broad 
|a(fc)| 2 while a long correlation length scale implies a narrow |a(/c)| 2 . Inserting Eq. ( |4.4| ) into 
Eq. ($LT§ yields 

A,(t) = \C I \a{k)\ 2 {l - \(e^)\ 2 )dkdr, (4.7) 



h jo 

where x{r) is given by Eq. ( |3.5|) , and e lkx ^ can immediately be recognized as Glauber's 
displacement operator, with a translation in momentum of hkcos(ur) and a translation in 
position of —{hk)/{mu) sin(c<jr). In terms of the translated state 

W;k,r) = e ikx V\4>), (4.8) 

the entropy production is given by 

M{t) = JiJoJ \< k )\^ l -M^ k ^)\ 2 )dkdT. (4.9) 

While it is not possible to find general solutions for the minimization of Eq. ( |4.9| ) for arbi- 
trary environment correlations, it is possible to extract important limiting cases. Restricting 
attention to times of several oscillator periods or more, entropy minimization requires the 
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maximization of \{ip\ip\k 1 r)\ 2 for typical values of k, on the order of 5k (the spread of 
|a(/c)| 2 ), resulting in typical translations of h5k in momentum and (h5k)/{muj) in position. 
The maximization of the square of the inner product of any two normalized vectors occurs 
when the vectors are identical (up to a phase). Thus the maximal states and the translated 
maximal states will be approximately equal for typical translations, requiring that the width 
of |a(A;)| 2 be much less than the width of the maximal states: 

5k < ——Ax 

h 

5k < ^Ap. (4.10) 

Since |a(A;)| 2 and c(r) are Fourier transform pairs, this last condition also implies that the 
noise spatial correlation function c(r) is wide (compared to the maximal states). In this 
long correlation length scale limit, g(x,x') in Eq. (|4.1|) is quadratic |pT[ , corresponding to 
the low dissipation limit studied by ZHP and to the results of Sec. |TT| with bj = for all j. 
The coherent states are then the maximal states if they are consistent with the condition 
expressed in Eq. ft4.10|) using coherent state values for Ax and Ap: 



which corresponds to a environment corellatation length much larger than the width of a 



coherent state Ax = yh/2muj. 

If the environment correlation length is shorter than the width of the coherent state, then 
the approximation described above is not valid. It is useful to examine the entropy produc- 
tion in the position representation and absorbing the time dependance into the (Schrodinger 
picture) state vectors, where 

{^\e ikx{r) \^) = (^\U\T)e ikx U{T)\^j) 

= mr)\e lkx mr)) 
dxe ikx \iP(x,r)\ 2 

dxe ikx P{x,T). (4.12) 
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Entropy production then becomes 

2 



Adt) = -n [c(0)t - / drdxdx 1 'c(x - x')P(x, t)P(x', t)}, (4.13) 

h J 



using the Fourier transform relation expressed in Eq. ( |4.6| ). For narrow c(r) this expression 
becomes independent of the function P. To see this, consider an example where the spatial 
correlation of the environment is given by 

c (r) = Ae { ^ )2 . (4.14) 

In the limit o approaches zero 

J dxdx'\e ( ^ )2 P(x,T)P(x',T) = 

XoyfH J P 2 (x,r)dx -> 0, (4.15) 

so that in the short correlation length limit 

Aq(t) = ^c(0)t (4.16) 

so that in this regime, all states produce the same entropy, there are no maximal states. 
However, the decay rate of the off diagonal terms given by the decoherence term in Eq. ( [Op 
in the low short correlation length regime is generally at a maximum, 

dp(x,x';t) 2c(0) 

= tf-p{z,x;t), (4.17) 

and if decoherence is to be effective, the decay time must be comparable to dynamical 
timescales. There will necessarily be a significant increase in the entropy for all pure states, 
so that all states will be rapidly ' mixed' by the noise. 

V. COMMENTS AND CONCLUSIONS 

I have established an approximation scheme for determining maximal states (as defined 
by ZHP), and applied it to two families of Lindblad operators. For Lindblad operators 
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which are at most quadratic in position and momentum, squeezed states emerge as the 
maximal states for intermediate times compared to the dynamical time scales. The amount 
of squeezing decreases with time, so that coherent states are maximal for large timescales. 
Large timescales are the most relevant, since an object's classicality should be an enduring 
property, not a transient one. For an environment with finite spatial correlation, coherent 
states emerge as maximal when the environment has long correlation length, but all states 
rapidly become mixed states when the environment length scale is long. Thus environment 
correlation effects will not be important in establishing the nature of maximal states and the 
character of quasiclassical states. However, correlation effects can still be important when 
considering quantum interference between two such states. One important result which 
emerges from these calculations is that coherent states are a robust choice for the maximal 
states. 
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